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Abstract. Observations of helium and hydrogen emission lines from metal-poor extragalac- 
tic H II regions, combined with estimates of metallicity, provide an independent method for 
determining the primordial helium abundance, Y p . Traditionally, the emission lines employed 
are in the visible wavelength range, and the number of suitable lines is limited. Furthermore, 
when using these lines, large systematic uncertainties in helium abundance determinations 
arise due to the degeneracy of physical parameters, such as temperature and density. Re¬ 
cently, Izotov, Thuan, & Guseva (2014) have pioneered adding the He I A10830 infrared 
emission line in helium abundance determinations. The strong electron density dependence 
of He I A10830 makes it ideal for better constraining density, potentially breaking the degen¬ 
eracy with temperature. We revisit our analysis of the dataset published by Izotov, Thuan, & 
Stasihska (2007) and incorporate the newly available observations of He I A10830 by scaling 
them using the observed-to-theoretical Paschen-gamma ratio. The solutions are better con¬ 
strained, in particular for electron density, temperature, and the neutral hydrogen fraction, 
improving the model fit to data, with the result that more spectra now pass screening for 
quality and reliability, in addition to a standard 95% confidence level cut. Furthermore, the 
addition of He I A10830 decreases the uncertainty on the helium abundance for all galaxies, 
with reductions in the uncertainty ranging from 10-80%. Overall, we find a reduction in the 
uncertainty on Y p by over 50%. From a regression to zero metallicity, we determine Y p = 
0.2449 ± 0.0040, consistent with the BBN result, Y p = 0.2470 ± 0.0002, based on the Planck 
determination of the baryon density. The dramatic improvement in the uncertainty from 
incorporating He I A10830 strongly supports the case for simultaneous (thus not requiring 
scaling) observations of visible and infrared helium emission line spectra. 
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1 Introduction 

Although the determination of the fundamental parameters of the standard cosmological 
model including dark matter and dark energy, ACDM, by WMAP [1, 2] and Planck [3, 4] are 

unparalleled, big bang nucleosynthesis and the observations of the light element abundances 
offer an important cross-check, in particular on the determination of the baryon density. The 
most recent Planck result [4] for the baryon density, H^/i 1 2 * = 0.02226 ± 0.00016, corresponds 
to a baryon-to-photon ratio of rj = (6.10 ±0.04) x 10 -10 . Because the uncertainty in r] is now 
less than 1%, standard big bang nucleosynthesis (SBBN) [5-7] is a parameter-free theory [8], 
and relatively precise predictions of the primordial abundances of the light elements D, 4 * * He, 

4 He, and 'Li are available [9-22], While the 7 Li abundance remains problematic [17], recent 
D/H determinations from quasar absorption systems have become quite precise, in their own 

right, and they are in excellent agreement with the prediction from SBBN and the CMB 
[23]. Using a neutron mean life of 880.3 ± 1.1s [24], SBBN yields a primordial abundance for 
4 He, Y p , of Y p = 0.2471 ± 0.0002, using the Planck determined value of 7] [22], By allowing 
Y p to vary as an independent parameter, fits to CMB anisotropies allow for a determination 
of Y p within the context of ACDM. The recent Planck results found Y p = 0.251 ± 0.027 [4], 
Fortunately, the helium abundance from emission line measurements provide significantly 
better precision. 

To test SBBN beyond D/H, it is clear that precise determinations of 4 He are necessary. 
In addition, 4 He still provides important constraints on the physics of the early universe 
beyond the standard model [25]. Nevertheless, obtaining better than 1% precision for in¬ 
dividual objects remains a challenge. 4 He abundance determinations are generally fraught 
with systematic uncertainties and degeneracies among the input parameters needed to model 
emission line fluxes [26, 27, 37]. The primordial abundance of 4 He is determined by fitting the 
helium abundance versus metallicity (oxygen), and extrapolating back to very low metallicity 
[28]- 

Our method for determining the 4 He abundance in an individual H II region is based 
on a Markov Chain Monte Carlo (MCMC) analysis [29-31]. Using calculated emissivities 
[32, 33], we model emission fluxes for a number of H and He lines relative to H/3. The model 
depends on several physical parameters associated with the H II region: electron density, 
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n e , temperature, T, optical depth r, underlying stellar H and He absorption, an and aue, 
reddening, C(H(3), the fraction of neutral hydrogen, and, of course, the He abundance, 
parametrized in terms of the abundance by number (relative to H) of ionized He, y + . MCMC 
scans of our 8-dimensional parameter space map out the likelihood function based on the y 2 
given by 

f FM _ EW \2 

2 _ ST Wmmeas’ n v 

x ^ o{\y ’ 1 j 

where the emission line fluxes, F( A), are measured and calculated for six helium lines (AA3889, 
4026, 4471, 5876, 6678, and 7065) and three hydrogen lines (Ha, Hy, H5), each relative to 
H/3. The x 2 i n eq. 1-1 runs over all He and H lines and u(A) is the measured uncertainty in 
the flux ratio at each wavelength. 

We also adopt a weak prior on the temperature based on the temperature derived 
from [O III] emission lines. Xt = (T ~ T(OIII)) 2 / a 2 ) is added to eq. 1.1, where we adopt 
(j = 0.2 T meas , which is much greater than the estimated temperature difference of 11% 
between the measurements of temperature from ionized helium and doubly ionized oxygen 
lines [34-36]. For more detail on the applied temperature prior see Aver et al. [29, AOS2]. 
Here we also discuss the implications of removing the prior. 

Once minimized, best-fit solutions for the eight physical parameter inputs, including 
y + , are found, and uncertainties in each quantity can be obtained by calculating a ID 
marginalized likelihood. With eight parameters and only nine observables, our model is 
nearly under-determined. 

This MCMC analysis was applied to the large HeBCD dataset of Izotov, Thuan, & 
Stasiriska [37, ITS07] in Aver et al. [30, AOS3]. As discussed in AOS3, we consider only 
those objects with measured He I A4026 emission lines. We then select those objects with 
X 2 < 4, corresponding to a standard ~95% confidence level for one degree of freedom. Using 
a revised set of emissivities [32, 33], in our most recent work, Aver et al. [31, AOPS], we 
found only 16 objects which satisfy this criterion and could be used to extrapolate to zero 
metallicity and determine Y p . Due to the large uncertainties in y + found for each object, and 
the relatively small number of objects for which the model acceptably describes the data, the 
uncertainty on our determination of Y p was not much better than 4%. 

Recently, Izotov, Thuan, & Guseva [38, ITG14] presented data for the near-infrared He I 
A10830 emission line for a large sample of metal-poor nebulae and pioneered the addition of 
this line in helium abundance analyses. While not only providing an additional test for the 
model, this particular line shows a strong sensitivity to the electron density. It can, therefore, 
be very useful in breaking the degeneracy between density and temperature [27, 39], allowing 
for more accurate determinations of y + in each region. Here, we make use of these new 
observations in a new MCMC analysis. We consider only observations for which all seven He 
lines have reported values in ITS07 & ITG14 [37, 38]. While the initial data set is reduced 
(not every object measured in ITS07 has a measured He I A10830 emission line in ITG14 and 
vice versa), better constrained parameters allow more objects to yield statistically reliable 
and physically meaningful solutions. As a result, we base our current analysis on 15 objects 
(with some, but not full, overlap with the sample used in AOPS [31]). The uncertainty 
in y + in individual objects is, in fact, significantly reduced (by as much as 80%) and the 
extrapolated primordial abundance is now determined to 1.6%: Y p = 0.2449 ± 0.0040. 

This paper is organized as follows. First, §2 introduces He I A10830, with particular em¬ 
phasis on its electron density dependence. Second, in §3, Monte Carlo testing with synthetic 
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data including He I A10830 is carried out to examine its diagnostic power. This includes ex¬ 
amining cases both with and without the inclusion of an 0 [III] temperature prior. Third, §4 
offers a brief overview of the ITG14 observations and dataset, along with the required scaling 
of the He I A10830 observations. In §5, the set of ITS07 galaxies in ITG14 are analyzed, 
demonstrating the impact of He I A10830 on the solutions. Subsequently, Y p is determined 
in §6. Finally, §7 offers a discussion of the results, possible systematic effects or additional 
uncertainties, and the benefits of further observations with He I A10830. 


2 The diagnostic power of the infrared He emission line, A10380 


He I A10830 is emitted in a 2 3 P to 2 3 S transition (A£ = 1). Corresponding to its small 
energy change (A E = 1.14 eV), it is readily collisionally excited from the metastable 2 3 S 
triplet state. It is this large population due to collisional excitation that makes it so sensitive 
to the H II region’s density. Figure 1 shows the sensitivity to density for all of the He emission 
lines considered. The figure clearly shows the dramatically stronger density dependence of 
He I A10830, with a slope over four times greater than that of the next most sensitive line, 
He I A7065. Therefore, we expect that He I A10830 should be an excellent density diagnostic, 
greatly constraining the density in our solutions. Figure 2 shows the temperature dependence 
of the He emission lines considered. He I A10830 does not stand out, as it did in figure 1, 
but its temperature dependence is the second strongest, and, furthermore, its increase with 
temperature helps to balance out the three strong optical lines, He I AA4471, 5876, & 6678, 
all of whose emissivities decrease similarly with temperature. 

The He line flux ratios (compared to H/3) used in eq. (1.1) are calculated using 


F( A) 
F(H/3) 


= y 


+ 


E{ A) 
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W (H fS )+ a H (H f 3 ) 
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1 O - f (\) C ( H 0 ) , 
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where W{ A) is the measured equivalent width, and two parameters, an and a# e , characterize 
the wavelength-dependent underlying absorption for H and He, respectively. The function 
f T ( A) represents a correction for radiative transfer and depends on on r, n e , and T. The 
emissivity, E, including the He collisional corrections, is a function of n e , and T. The 
hydrogen collisional corrections are accounted for by which depends on the fraction of 
neutral-to-ionized hydrogen, £. The final term in eq. (2.1) accounts for reddening, C(H/3). 
See AOS &; AOS2 [29, 39] for full details on the model and analysis method. 

To incorporate He I A10830 into our analysis, the following data are used. The emissivi¬ 
ties, including the collisional correction, are adopted from the recent work of Porter, Ferland, 
Storey, & Detisch [33], as was implemented in AOPS [31], and include He I A10830. The 
radiative transfer equations in terms of optical depth come from the work of Benjamin, Skill- 
man, & Smits [40]. That work does not include a fitting formula for He I A10830, however. 
Instead, the 10-level numerical calculation program was graciously provided by Bob Ben¬ 
jamin (private communication). Fitting the same functional form to the He I A10830 optical 
depth data over the same range (n e = 1-300 cm~ 3 and T = 1.2-2.0 x 10 4 K), we obtain: 


/ r ( 10830) = 1 + (t/ 2)[0.0149 + (4.45 x 10“ 3 - 6.34 x 10" 5 n e + 9.20 x 10” 8 n e 2 )T)], (2.2) 


with a maximum fit error of 1.3% (similar to the maximum fit errors reported in [40]). 
Correcting for the effects of underlying absorption for He I A10830 is more tenuous since the 
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Figure 1. Comparison of the Porter, Ferland, Storey, & Detisch [33] emissivities, including the 
collisional correction, versus density, using a temperature of 18,000 K, for He I AA3889, 4026, 4471, 
5876, 6678, 7065, 10830. The much stronger density dependence of He I A10830 is apparent. 


available data is very limited. The 3 B-star supergiant spectra (where helium absorption is 
more significant) reported in Conti & Howarth [41] show equivalent width (EW) underlying 
absorption of ~0.5-2.0 A. Because the EW of He I A10830 is large (typically 100-400 A), 
the correction of underlying absorption will be reassuringly small. For this work, we adopt 
an underlying absorption coefficient of 0.8, relative to He I A4471 (i.e., He I A4471 carries a 
coefficient of 1.0), such that the underlying absorption correction applied will be 0.8 x a# e , 
where aue is the solution for underlying absorption (see AOS for further details [39])). This 
ratio of 0.8 for underlying absorption of He I A10830 relative to that of He I A4471 is the 
same as in adopted by ITG14 [38] (though the approaches to determining the underlying 
absorption differ). 
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Figure 2. Comparison of the Porter, Ferland, Storey, & Detisch [33] emissivities, including the 
collisional correction, versus temperature, using a density of 100 cm -3 , for He I AA3889, 4026, 4471, 
5876, 6678, 7065, 10830. 

3 Characterizing the impact of He I A10380 with synthetic testing 

To begin discerning the effect and potential benefits of He I A10830, we conduct Monte Carlo 
analyses using synthetic data, since it offers a controlled environment for clear comparisons. 
First, a set of input parameters were used to generate a synthetic spectrum. These are given 
in the second column of table 1. Then, 1000 sets of Gaussian perturbed fluxes were calculated 
with spreads of 3% and 1.5% of the unperturbed helium and hydrogen fluxes, respectively. 
The best-fit solution (via y 2 minimization, see AOS [39]) of each of these 1000 perturbed 
spectra was then found. Finally, the average value and dispersion of those 1000 solutions was 
calculated for each parameter. To allow comparisons showing the effects of He I A10830 and 
the [O III] temperature prior, this process was repeated for the following four cases: analysis 
with He I A10830 & with T(0 III) included, without He I A10830 but with T(0 III), with 
He I A10830 but without T(0 III), and without He I A10830 & without T(0 III). The results 
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of these MC simulations for y + vs. both n e and T (the electron density and temperature 
solved for from the helium and hydrogen emission lines) are shown in figure 3. The average 
value and dispersion for each of the four cases is provided in table 1. 




w/ T(0 III) & 



w/o T(0 III) & 


Input 

w/ He I A10830 

w/ T(0 III) 

w/ He I A10830 

w/o He I A10830 

y + 

0.085 

0.0858 ± 0.0027 

0.0857 ± 0.0063 
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a He 

0.5 

0.52 ± 0.11 

0.51 ± 0.15 

0.51 ± 0.12 

0.53 ± 0.14 

T 

1.0 

0.78 ± 0.31 

0.80 ± 0.73 

0.62 ± 0.56 

0.52 ± 0.72 

T e 

16,000 

17,320 ± 1090 

17,400 ± 1300 

18,860 ± 4190 

18,300 ± 5320 

C(H/3) 

0.1 

0.09 ± 0.03 

0.09 ± 0.03 

0.09 ± 0.03 

0.08 ± 0.04 

a h 

1.0 

2.08 ± 1.47 

2.10 ± 1.70 

2.61 ± 2.09 

2.69 ± 2.15 

i x io 4 

1.0 

i o +19 
—13 

14 +31 
_ 14 

1040 +31,62 ° 

-Ll+lU -1040 

ocyn +68530 

OO/U _ 85 7o 

T(0 III) 

17,000 






Table 1. Comparison of MC distributions with synthetic data 


When the [O III] prior and He I A10830 are included, the input data are nicely repro¬ 
duced. As the figure and table clearly show, He I A10830 dramatically reduces the uncertainty 
on the density (by over 75%), as well as improving the average, both with and without the 
inclusion of the T(0 III) prior. As a result of the degeneracy between temperature and 
density (see OS04 &; AOS [27, 39]), the uncertainty on the temperature is also decreased. 
Corresponding to the better constrained density and temperature, the precision on the helium 
abundance, y + , is greatly increased. The inclusion of He I A10830 decreases the uncertainty 
on y + by over 50%, when either the T(0 III) prior is or is not included. In fact, in terms of the 
uncertainty on y + , the inclusion of He I A10830 mostly obviates the benefit of the T(0 III) 
prior. However, the T(0 III) prior still has a large effect in improving the temperature 
determination. 

The neutral hydrogen fraction, £, is very poorly constrained for very low temperatures, 
due to its exponential temperature dependence (see AOS [39]). As a result, the T(0 III) 
prior’s exclusion of the lowest temperatures is crucial in restraining the range of MC solu¬ 
tions for £. Objects with completely unphysical best-fit solutions for £ (i.e., >25% neutral 
hydrogen) are excluded (see AOS3 [30]). Therefore, the T(0 III) prior still offers potential 
utility in our analysis here, even if it does not improve the determination of y + when He I 
A10830 is included. It should be noted that we first implemented T(0 III) as a weak prior 
to rule out unphysical double minima in some solutions (see AOS2 [29]). The much better 
constrained density in the presence of He I A10830 has a similar effect, thus allowing for the 
possibility of removing the non-parametric prior. 

As can be seen particularly in the lowest of the density panels in figure 3, some MC 
realizations have solutions for which the density increases to very large values. This is a con¬ 
sequence of the temperature-density degeneracy, and as discussed above, the corresponding 
very low temperatures allow the neutral hydrogen fraction to increase almost without bound 
(see AOS for further discussion [39]). To see the extent of the range in density found in this 
case, we provide an inset in that figure. As a result of the extended range in n e , the MC 
average values are biased upward for n e and £, with a correspondingly large uncertainty, as 
can be seen in the inset and the last column of table 1. This is one of the drawbacks of Monte 
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Figure 3. Comparison of 1000 MC solutions for analysis with and without He I A10830, as well as 
with and without the [O III] temperature prior. For each of these four cases, the helium abundance 
is plotted separately versus density and temperature. On each of these eight panels the average MC 
best-fit value is plotted with its dispersion. From top to bottom the progression of the panels is as 
labeled: w/ He I A10830 & w/ T(0 III), w/o He I A10830 & w/ T(0 III), w/ He I A10830 & w/o 
T(0 III), and w/o He I A10830 & w/o T(0 III). The panels on the left are y + vs. n e , while those on 
the right are y + vs. T, and, except for the inset in the bottom left panel, all of the axes are shared 
and show the same domain and range. 
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Carlo via flux perturbation, as was discussed in AOS [39]. The MCMC method introduced in 
AOS2 [29] is not sensitive to this same behavior. Nevertheless, the distribution of perturbed 
solutions is very effective at illustrating the effects of the T(0 III) prior and of He I A10830. 

For a further demonstration of the effects of He I A10830, one of the flux perturbed 
synthetic spectra was analyzed using the MCMC analysis we employ to determine the so¬ 
lutions to the physical parameter set [29-31]. Rather than examining the distribution of 
1000 flux simulations, the spectrum, synthetic or otherwise, is treated as an observation, for 
which a Monte Carlo over the eight model parameters (y + , n e , a# e , r, T, C(H/3), a h, £) 
is performed to determine the uncertainties about the best-fit values. For one of the flux 
perturbed synthetic spectra, the solution (with uncertainties) was found for the same four 
cases discussed above and is reported in table 2. The conclusions from the above analysis 
are reinforced, with He I A10830 tremendously improving the density determination, both 
the best-fit value and the uncertainty. This, in turn, improves the helium abundance de¬ 
termination. Comparing the solutions with and without He I A10830, the 68 % confidence 
level uncertainty range for the density is decreased by ~75% &; ~80%, for the cases with 
and without the T(0 III) prior, respectively, while the uncertainty in y + sees corresponding 
decreases of ~40% &; ~50%. Again, the efficacy of He I A10830 in improving our solutions 
is impressive. 
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3.11 
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Table 2. Comparison of MCMC solutions with synthetic data 


4 The ITG14 observations and dataset 

Izotov, Thuan, & Guseva [38, ITG14] have recently published their observations of helium 
abundances using He I A10830. The ITG14 dataset consists of near-infrared measurements of 
He I A10830 and hydrogen emission line Paschen 7 A10940 A (P 7 ) from 45 H II regions in 43 
low metallicity galaxies. The majority of the observations, 33 of 45, were made with the 3.5 
m Apache Point Observatory and TripleSpec spectrograph; 8 observations were taken with 
the 8.4 m Large Binocular Telescope and Lucifer spectrograph, as it allows observations of 
fainter, often lower metallicity, galaxies; and finally, four observations were retrieved from the 
European South Observatory data archives (2 VLT & 2 NTT). ITG14 combine these near- 
infrared spectra with 75 optical spectra taken from those 43 galaxies, where some galaxies 
have multiple independent observations. Because the near-infrared and optical spectra were 
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taken with different telescopes, apertures, and at different times, the measured He I A10830 
flux must be scaled to match the optical flux measurements. Similar to ITG14 [38], we scale 
He I A10830 using the observed P 7 flux and the theoretical flux ratio for P 7 relative to H/3. 


F(A10830) 

F(Hp) 

F( A10830) 
F(H0 ) 


F(A10830) F(P'y) 
F(P^) F(H{3) 


^ W{Hp)+a H {Hp) 

\10830) E(P 7 ) w(Hfi) 


(4.1) 


The theoretical flux ratio for P 7 relative to H/3 is calculated analogously to eq. (2.1). 
As for all of the hydrogen emission lines, the emissivity ratio for P 7 relative to H/3 is taken 
from Hummer & Storey [42]. The collisional correction for P 7 is calculated, as in AOS [39], 
using collisional excitation rates from Anderson et al. [43] (similar to H<3) and recombination 
rates from Hummer & Storey [42], The estimate of the underlying absorption correction for 
P 7 is potentially much more significant than for He I A10830 because it is an intrinsically 
weaker emission line and is constrained by the same lack of data as in the case of He I A10830. 
The same 3 B-star supergiants discussed in §2 show a similar ~0.5-2.0 A EW of underlying 
absorption for P 7 [41]. For this work, we adopt an underlying absorption coefficient of 0.4, 
relative to H/3 (i.e., H/3 carries a coefficient of 1.0), such that the underlying absorption 
correction applied will be 0.4 x an, where an is the solution for underlying absorption based 
on all of the hydrogen lines (see AOS for further details [39]). 

Of the 45 H II regions with measured He I A10830 emission lines reported in ITG14 
[38], 26 are also present in the 93 observations comprising the HeBCD dataset of ITS07 [37], 
which we have previously analyzed in AOS3 &; AOPS [30, 31]. Of the 93 HeBCD objects, 
He I A4026 is detected for 70. To reduce systematic uncertainty that may be introduced by 
the absence of He I A4026 (see AOS3 for further details [30]), objects where He I A4026 is not 
reported are excluded in our analysis. Of the 26 H II regions in the HeBCD for which He I 
A10830 is reported in ITG14, He I A4026 is detected for 22. However, of those 22, several 
have multiple independent observations, increasing the total number of spectra with both 
He I A10830 and He I A4026 to 31. Those 31 sets of observations, with optical spectra taken 
from ITS07 and near-infrared from ITG14, comprise our initial dataset for this work. 


5 Tracking the effects of He I A10830 

For the 31 HeBCD observations from ITS07 [37] with He I A10830 observations reported 
in ITG14 [38], the effects of including He I A10830 into our helium abundance analysis is 
examined for two cases. First, we examine the effects of the addition of He I A10830 to the 
analysis of AOPS [31]. As discussed earlier, in AOS 2 , AOS3, & AOPS [29-31], a temperature 
prior based on [O III] emission lines was utilized to eliminate unphysical double minima and 
improve the parameter determination. However, it is a non-parametric (i.e., not fully “self- 
consistent”) portion of our model, and it introduces a small bias in the solution. Due the 
inclusion of He I A10830, the T(0 III) prior is no longer needed to eliminate double minima, 
and synthetic testing in §3 showed that its benefit in reducing parameter uncertainties was 
greatly reduced by the inclusion of He I A10830. Therefore, we also examine the effects of 
removing the T(0 III) prior once He I A10830 is included. 
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For each of these analyses, cuts on the dataset are made following AOS3 [30] (also 
employed in AOPS [31]). Of primary interest is the standard 95% confidence level cut. In 
AOS3, there was one degree of freedom (9 fluxes relative to H (5 and 8 model parameters). 
However, with the addition of He I A10830, there are now two degrees of freedom in our 
analysis. Therefore, a 95% confidence level corresponds to x 2 < 6 , and best-fit solutions with 
X 2 > 6 were excluded. Furthermore, objects with unphysical physical parameters, specifically 
£ > 0.333 (>25% neutral hydrogen), were also excluded. Of the remaining, qualifying objects, 
objects with large corrections for systematic effects were flagged to limit potential systematic 
errors introduced by uncertainties in the models. Objects with r > 4 and £ > 0.01, where 
the 1-a lower bound does not encompass £ = 0 . 001 , were flagged. 

In the analysis of AOPS [31], the 31 objects in our present dataset broke down as follows: 
11 qualifying, 2 flagged (for r > 4), and 18 excluded, 12 for x 2 > 4, 3 for £ > 0.333, and 
3 for both. When He I A10830 is added, the 11 previously qualifying objects are retained 
as qualifying. One of the flagged objects is now qualifying, the other is still flagged for r. 
Six previously excluded objects are no longer excluded, four qualifying and two flagged (for 
£). Thus, the number of qualifying objects increases from 11 to 16, while the total number 
of flagged objects increases by one (from 2 to 3). If the T(0 III) prior is then removed, 
only three originally excluded objects are reclaimed, one qualifying, and two flagged (for £), 
and one originally flagged object is excluded (the other remains flagged for r). As a result, 
the qualifying and flagged datasets each increase by only one, from 11 to 12 and 2 to 3, 
respectively. As a final note on qualifying and flagged objects, two of the qualifying and 
one of the flagged observations in the analysis retaining the T(0 III) prior, SBS 0335-052E1, 
SBS 0335-052E3, & SBS 0335-052E3, respectively, are independent observations of the same 
H II region. As such, they are combined by weighted average into a single object, reducing 
the number of qualifying objects from 16 to 15, and the distinct flagged objects from 3 to 2. 

For the 11 jointly qualifying objects in AOPS and this work, figures 4 & 5 compare the 
best-fit solutions and uncertainties for y + , n e , and T 4 (T e /10 4 ), first for the addition of He I 
A10830, and then for the removal of the T(0 III) prior. In figure 4, the effect of including 
He I A10830 is examined. The AOPS values (whose analysis utilized the T(0 III) prior) 
are plotted adjacent to the results where He I A10830 is included (and the T(0 III) prior is 
retained). The uncertainties for y + show decreases ranging from 10-80% (see upper panel). 
As was seen in §3, this substantial improvement in the determination of y + stems from the 
dramatically better constrained density. The 68 % confidence level density range is reduced 
by over 60%, with most objects seeing their uncertainty range on the density decrease by 
around 85% (see middle panel). Due to the temperature-density degeneracy (see OS04 & AOS 
[27, 39]), the better constrained density results in a reduced temperature uncertainty (see 
lower panel). As density and temperature are the most important parameters in relating the 
helium abundance to the observed flux, y + is correspondingly better determined. Note that 
in all cases, the new solutions are fully consistent with the prior solutions, though now with 
greatly reduced uncertainties. The best-fit values for y + decrease overall (4 of 11 increase, 7 
decrease), but the average reduction in the helium abundance is small, approximately 2 %. 

Figure 5 shows the same parameter comparisons but examines the effect of removing 
T(0 III) prior after He I A10830 has been added. Therefore, both series include He I A10830, 
and the results without and with the T(0 III) are plotted adjacent to each other. Overall, 
there is a slight improvement in the parameter determinations from the inclusion of the 
T(0 III) prior, but it is not nearly as dramatic as in figure 4. The uncertainties on y + 
generally decrease by 5-10% due to the inclusion of the T(0 III) prior (see figure 5, upper 
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panel). 

Because the addition of He I A10830 increases the constraint on the density so strongly, 
the T(0 III) prior is no longer needed to eliminate unphysical double minima, while the 
T(0 III) prior’s benefit in determining y + is greatly reduced (as was shown in §3 and figure 
5). However, as discussed in §3, the neutral hydrogen fraction, £, is very poorly constrained 
for low temperatures, due to its exponential temperature dependence (see AOS [39]). As a 
result, the T(0 III) prior combined with He I A10830 helps to constrain the temperatures 
further than either does individually. The preclusion of the lowest temperatures restricts £ 
and prevents the solution being excluded due to a completely unphysical best-fit value for £ 
(>25% neutral hydrogen). Therefore, because the T(0 III) prior allows three more objects 
to qualify, increasing the qualifying dataset by 25% compared to adding He I A10830 but 
removing the T(0 III) prior, we favor retaining the T(0 III) prior for this analysis. 

This dataset of 15 qualifying objects and 2 flagged objects (analyzed including both 
He I A10830 &; the T(0 III) prior) comprises our Final Dataset. These are the objects for 
which the model is a good fit and which returns physically meaningful parameter solutions. 
As discussed above and as shown in figure 4, the parameter determinations of their solutions 
are significantly improved compared to our analysis in AOPS [31], which did not include He I 
A10830. The best-fit solutions and uncertainties of the Final Dataset are presented in table 
3, and they are used to determine Y p in the following section (§6). 
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Figure 4. Plot comparing the parameter solutions for y + , n e , and X 4 (T e /10 4 ) for the 11 qualifying 
objects in AOPS with He I A10830 values in ITG14. The lighter, thinner lines show the results 
given in AOPS. The results after the inclusion of He I A10830 are given in the darker, thicker bars. 
The uncertainty in n e is dramatically reduced by the addition of He I A10830 with a corresponding 
reduction in the uncertainty in y + . 
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Figure 5. Plot comparing the parameter solutions for y + , n e , and X 4 (T e /10 4 ) for the 11 qualifying 
objects in AOPS with He I A10830 values in ITG14. The lighter, thinner lines show the results where 
He I A10830 is included but the T(0 III) prior is removed. The results including both He I A10830 
and the T(0 III) prior are given in the darker, thicker bars. Because including the T(0 III) prior does 
not significantly bias the values of y + , and because its inclusion results in a greater yield of qualifying 
points, the T(0 III) prior is applied in the final analysis. 
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Table 3. Physical parameters and He + /H + abundance solutions of the Final Dataset 









6 Results from the Final Dataset 


Following the results of the previous sections, we now calculate the primordial helium abun¬ 
dance (mass fraction), Y p , from the data listed in table 3. A regression of Y, the helium mass 
fraction, versus O/H, the oxygen to hydrogen mass fraction, is used to extrapolate to the 
primordial value 1 . The O/H values are taken directly from ITS07[37]. 

The relevant values for the Final Dataset’s regression are given in table 4. The regression 
based on the 15 qualifying objects yields, 

Y p = 0.2449 ± 0.0040, (6.1) 

with a slope of 79 ± 43 and a total y 2 of 7.6. The result is shown in figure 6 and summarized 
in table 5. This result for Y p agrees well with the SBBN value of Y p = 0.2471 ± 0.0002 [22], 
based on the Planck determined baryon density [4], Not surprisingly, eq. (6.1) also agrees 
well with the SBBN-independent, direct Planck estimation of Y p = 0.251 ± 0.027. AOPS 
determined Y p = 0.2465 ± 0.0097 with a slope of 96 ± 122. The result from AOPS is in good 
agreement with the newer result. Restricting the qualifying AOPS dataset to the 11 objects 
which are shared with the qualifying dataset of this work, results in Y p = 0.2461 ± 0.0107 
with slope 84 ± 144, again in good agreement with this result found in eq. (6.1). Due to five 
more objects qualifying, the Final Dataset of this work contains a similar number of objects 
as in AOPS (16 there, 15 here). However, as a result of the reduction in the individual 
y + uncertainties (see §5), the uncertainty on the intercept and slope are each reduced by 
approximately 60%. The intercept value decreased by 0.65% compared to AOPS. Given 
that the uncertainty on the intercept determination in AOPS was 3.9%, this decrease is not 
significant. 

Including the 3 flagged objects 2 decreases the intercept and reduces the uncertainty 
to 0.2424 ± 0.0034 with a slope of 116 ± 32. The reduced uncertainty is primarily a re¬ 
sult of the increased number of points in the regression, while the decreased intercept is 
entirely the result of the increased slope. As was seen in AOS3 &; AOPS, the flagged data 
points tend to have higher helium abundances, primarily due to the objects flagged for large 
neutral hydrogen fractions. In previous analyses, where the evidence for a non-zero slope 
was marginal, we have reported the mean value of Y. The mean value of Y for the Final 
Dataset using only qualifying points is <Y>= 0.2515 ± 0.0017. The mean value increases 
to <Y>= 0.2533 ± 0.0016, when the flagged objects are included. Olive & Skillman [27] re¬ 
stricted the metallicity baseline to O/H < 9.2 x 10~ 5 . Adopting the same metallicity cut with 
the dataset of this work increases the intercept to 0.2466 ± 0.0063, in near perfect agreement 
with the SBBN result using the Planck baryon density. Note, however, that in this case the 
reduced baseline leaves us with an undetermined slope, and for this case alone, basing the 
primordial abundance on the mean value given in the table 5 is justified. 

Instead, if a regression analysis is performed excluding the T(0 III) prior, the 12 qual¬ 
ifying objects yield Y p = 0.2379 ± 0.0050. However, the loss of the information provided 
by the T(0 III) prior and the reduced sample size weakens the value of this determination. 
Nevertheless, the substantial reduction in the value of the intercept raises some questions 
regarding possible systematic errors when the prior is dropped. 

1 This work takes Z = 20 (O/H) such that Y = 

2 One of the flagged objects is another independent observation of SBS 0335-052E that is combined with 

its two qualifying observations. 
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Figure 6. Helium abundance (mass fraction) versus oxygen to hydrogen ratio regression calculating 
the primordial helium abundance. 


Starting with the 75 observations corresponding to the 45 H II regions for which they ob¬ 
served near-infrared spectra to measure He I A10830, ITG14 select 28 for which W ( Hj3 ) > 150 
A, O++/0 > 0.8, and a(Y)/Y < 3%. Their resulting regression analysis and accounting of 
the dominant systematic uncertainties finds Y p = 0.2551 ± 0.0022. The significant differences 
between their result and the result derived here are due both to differences in determining 
the helium abundances for individual objects and in how the samples are defined. 

7 Discussion 

The primary aim of this work was to evaluate the effects of adding He I A10830 to a helium 
abundance analysis based solely on optical emission lines. In particular, it was hoped that 
this emission line, due to its much stronger electron density dependence, would have a sig¬ 
nificant effect in constraining the density in our solutions and, thereby, effectively break the 
degeneracy between temperature and density. Preliminary synthetic testing and analysis of 
H II region observations showed the positive impact of He I A10830 in improving the deter¬ 
mination of our solution parameters, n e , T e , and of primary interest, y + . Most objects saw 
their y + uncertainty decrease by ~50%, with a corresponding, similar improvement in the 
precision of Y p . 

Our determination of Y p = 0.2449 ± 0.0040 is in good agreement with the SBBN pre¬ 
dicted value (as well as our previous results), and while the value of the intercept has de¬ 
creased, it differs from the SBBN prediction by approximately 1/2 a, even with its most wel- 
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Object Ho + /H^ Hc++/H+ Y O/H x W 


Final Dataset Not Flagged (Qualifying) 


I Zw 18 SE 1 

0.07693 

± 

0.00423 

0.0008 

± 

0.0008 

0.2371 

± 

0.0100 

1.5 ± 0.1 

SBS 0335-052El+3< a ) 

0.08201 

± 

0.00303 

0.0026 

± 

0.0018 

0.2524 

± 

0.0080 

2 ± 0.1 

J0519+0007 

0.08875 

± 

0.00461 

0.0021 

± 

0.0021 

0.2665 

± 

0.0109 

2.8 ± 0.1 

SBS 0940+544 2 

0.08179 

± 

0.00316 

0.0005 

± 

0.0005 

0.2476 

± 

0.0073 

3.2 ± 0.1 

Tol 65 

0.07883 

± 

0.00294 

0.0011 

± 

0.0011 

0.2421 

± 

0.0072 

3.5 ± 0.1 

SBS 1415+437 (No. 1) 3 

0.07694 

± 

0.00494 

0.0022 

± 

0.0022 

0.2402 

± 

0.0125 

4 ± 0.1 

SBS 1415+437 (No. 2) 

0.08226 

± 

0.00343 

0.0000 

± 

0.0000 

0.2474 

± 

0.0078 

4.2 ± 0.3 

CGCG 007-025 (No. 2) 

0.08867 

± 

0.00623 

0.0007 

± 

0.0007 

0.2631 

± 

0.0136 

5.5 ± 0.2 

Mrk 209 

0.08114 

± 

0.00272 

0.0011 

± 

0.0011 

0.2471 

± 

0.0066 

6.1 ± 0.1 

SBS 1030+583 

0.07937 

± 

0.00306 

0.0021 

± 

0.0021 

0.2454 

± 

0.0084 

6.4 ± 0.2 

Mrk 71 (No. 1) 

0.08539 

± 

0.00445 

0.0008 

± 

0.0008 

0.2560 

± 

0.0100 

7.2 ± 0.2 

SBS 1152+579 

0.08139 

± 

0.00303 

0.0012 

± 

0.0012 

0.2478 

± 

0.0073 

7.7 ± 0.2 

Mrk 59 

0.08548 

± 

0.00405 

0.0010 

± 

0.0010 

0.2564 

± 

0.0092 

10.1 ± 0.2 

SBS 1135+581 

0.08462 

± 

0.00072 

0.0008 

± 

0.0008 

0.2542 

± 

0.0025 

11.7 ± 0.3 

Mrk 450 (No. 1) 

0.08634 

± 

0.00441 

0.0003 

± 

0.0003 

0.2565 

± 

0.0097 

15.2 ± 0.4 


Final Dataset with Flags 

SBS 0335-052El+2+3^ a l 
Mrk 1315 

Mrk 1329 

0.08132 ± 0.00243 0.0025 ± 0.0015 

0.09341 ± 0.00274 0.0000 ± 0.0000 

0.09060 ± 0.00425 0.0000 ± 0.0000 

0.2509 ± 0.0064 
0.2710 ± 0.0058 
0.2650 ± 0.0091 

2 ± 0.1 

18.9 ± 0.4 
19.2 ± 0.5 


( a ) SBS 0335-052E1, SBS 0335-052E2, & SBS 0335-052E3 are independent observations of the same 
H II region. As such, they were combined by weighted average into a single regression point. 


Table 4. Primordial helium regression values 


Dataset 

N 

Y p 

dY/d(0/H) 

<Y> 

Qualifying 

15 

0.2449 ± 0.0040 

79 ± 43 

0.2515 ± 0.0017 

Qualifying + Flagged 

17 

0.2424 ± 0.0034 

116 ± 32 

0.2533 ± 0.0016 

Qualifying^), O/H < 9.2 x 10" 5 

12 

0.2466 ± 0.0063 

35 ± 125 

0.2482 ± 0.0025 

Qualifying, T(0 III) prior removed 

12 

0.2379 ± 0.0050 

145 ± 52 

0.2508 ± 0.0019 


( a ) Olive & Skillman [27] restricted the metallicity baseline to O/H < 9.2 x 10 5 , and the same 
dataset was analyzed in AOS & AOS2 [29, 39]. 


Table 5. Comparison of Y p for selected datasets 


come increased precision. There are, however, potential systematic effects to be addressed. 
First, the inclusion of the near-infrared He I A10830 with existing optical spectra requires 
scaling of He I A10830 (see §4). This naturally raises the possibility of introducing additional 
error or systematic bias. Second, as discussed in Porter et al. [44] and updated in Porter, 
Ferland, Storey, &; Detisch [32], the helium emissivities, the primary driver in determining 
the helium abundance from the observed fluxes, carry estimated uncertainties of ~0.2-0.7%. 
Furthermore, significant, systematic improvements in the helium emissivities have been made 
in recent years [32, 45]. Any further shifts in these emissivities will have a significant effect 
on the determination of Y p . 

The impressive benefits from including He I A10830 in nebular helium abundance de¬ 
terminations strongly supports the case for additional near-infrared spectral observations of 
low-metallicity H II regions. In particular, to remove potential systematic effects introduced 
by scaling He I A10830 observations, simultaneous observations of the optical and infrared 
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spectra would be decisive. Coupled with the opportunity for higher quality, higher resolu¬ 
tion spectra, collection of simultaneous optical and infrared spectra would be of considerable 
value and utility. In sum, there is potential to significantly reduce the uncertainty in the 
determination of Y p . 
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